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Abstract 

Wc report on the discovery of a diffuse X-ray source with ASCA, presumably associated with a molecular 
cloud in the vicinity of the supernova remnant RX J1713.7— 3946. The energy spectrum (1-10 keV) of 
the hard X-ray source shows a flat continuum, which is described by a power-law with a photon index of 
T = l.Olg'g. We argue that this unusually flat spectrum can be best interpreted in terms of characteristic 
brcmsstrahlung emission from the loss-flattened distribution of cither sub-relativistic protons or mildly 
relativistic electrons. The strong shock of RX J1713.7— 3946, which is likely to interact with the molecular 
cloud, as evidenced by CO-line observations, seems to be a natural site of acceleration of such nonthermal 
particles. The observed luminosity of Lx = 1.7 x 10 35 erg s" 1 (for a distance of 6 kpc) seems to require 
a huge kinetic energy of about 10 50 erg in the form of nonthermal particles to illuminate the cloud. The 
shock-acceleration at RX J1713.7— 3946 can barely satisfy this energetic requirement, unless (i) the source is 
located much closer than the preferred distance of 6 kpc and/or (ii) the mechanical energy of the supernova 
explosion essentially exceeds 10 51 erg. Another possibility would be that an essential part of the lost energy 
due to the ionization and heating of gas, is somehow converted to plasma waves, which return this energy 
to nonthermal particles through their turbulent reacceleration on the plasma waves. 
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ISM: supernova remnants — X-rays: individual (RX J1713.7— 3946) 



1. Introduction 

Supernova remnants (SNRs) are commonly believed to 
be major sites for the production of Galactic cosmic rays. 
The diffusive shock acceleration mechanism naturally ac- 
counts for the hard power-law production spectra of cos- 
mic rays in their sources, but, as long as the particle 
injection rate remains an unsolved problem, the theory 
docs not tell us conclusively what fraction of the initial 
kinetic energy of an explosion can be transferred to cosmic 
rays (see e.g. the recent review by Malkov, Drury 2001). 
Therefore, it is extremely important to derive the total en- 
ergy contained in subrelativistic and relativistic particles 
from observations of the relevant components of nonther- 
mal electromagnetic radiation. 

Radio synchrotron emission provides information about 
GeV electrons in most shell-type SNRs. Synchrotron X- 
ray emission recently discovered in SN 1006 (Koyama et 
al. 1995) and some other shell-type SNRs indicates that 
the electron acceleration continues effectively up to multi- 
TcV energies. It is believed, however, that the acceler- 
ated protons constitute the major fraction of nonthermal 
energy to which the mechanical energy of a supernova 
explosion is transferred through shock acceleration. The 
best way to explore the proton component of accelerated 



particles in SNRs is the detection of gamma-rays by the 
decay of ir° mesons produced in collisions between cosmic- 
ray protons and ambient matter (Drury, et al. 1994; Naito, 
Takahara 1994). The detection of gamma-rays of hadronic 
origin from SNRs by current and planned space- and 
ground-based instruments requires high-density environ- 
ments in or close to the particle-acceleration sites. Large 
molecular clouds interacting with shells of SNRs may act 
as effective gas targets for the production of tt° mesons and 
their subsequent decay to high-energy 7-rays (Aharonian, 
et al. 1994). It has recently been claimed (Butt et al. 2001) 
that there exists a compelling positional association of 
the unidentified 7-ray source 3EG J1714— 3857 with a 
dense molecular cloud (cloud A) overtaken, most prob- 
ably, by the shock front of the SNR RX J1713. 7-3946 (G 
347.3—0.5). Moreover, it has been argued that the TeV 
gamma-rays from this SNR detected by the CANGAROO 
collaboration are also the result of interactions of accel- 
erated protons with nearby dense gas targets (Enomoto 
et al. 2002), though Reimer and Pohl (2002) claimed that 
this interpretation seems to be inconsistent with EGRET 
observations of this region. 

The flux of the subrelativistic component of the accel- 
erating protons is generally inconclusive both on obser- 
vational and theoretical grounds, but it is quite possible 
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that this component dominates the total energetics of non- 
thermal particles. Unfortunately, subrelativistic protons 
arriving at the Earth are not a representative sample of 
the cosmic-ray flux in the Milky Way owing to the effect 
of solar modulation. Ionization losses of these particles 
play an important role in heating both diffuse neutral gas 
and molecular clouds, and in generating free electrons in 
molecular clouds, which are crucial for interstellar chem- 
istry. It is therefore of great importance to measure the 
energy released in the form of subrelativistic particles by 
cosmic-ray accelerators. 

The ideal diagnostic tool to probe subrelativistic cos- 
mic rays is the nuclear prompt 7-ray line emission at McV 
energies (Ramaty, Lingenfelter 1979). However, conserva- 
tive estimates show that the fluxes of even most promi- 
nent 7-ray lines are only marginally detectable because 
of limited sensitivities of gamma-ray instruments in the 
MeV band. On the other hand, as we discuss below, 
the search for the sites of concentrated subrelativistic pro- 
tons using bremsstrahlung X-rays could be very promis- 
ing with the help of superior sensitivities of current X- 
ray detectors, even though only a small fraction (~ 10~ 5 ) 
of the nonthermal energy of subrelativistic protons and 
electrons is radiated away in the form of bremsstrahlung 
X-rays. Remarkably, if the subrelativistic cosmic rays 
reside in a dense gas environment, the ionization losses 
would result in a significant flattening of the low-energy 
spectra, and thus leading to the characteristic "1/e" type 
bremsstrahlung X-ray spectrum (Uchiyama et al. 2002). 
Such an unusually hard X-ray spectrum may serve as a 
distinct signature of nonthermal bremsstrahlung origin of 
X-ray emission. 

In our previous paper (Uchiyama et al. 2002) we re- 
ported on the detection of hard X-rays from a localized 
region in the SNR 7 Cygni, which we attributed to 1/e 
bremsstrahlung from the loss-flattened electron distribu- 
tion. In this paper we present a more pronounced example 
of such hard X-radiation arriving from a massive cloud in 
the vicinity of SNR RX J1713. 7-3946. 

2. SNR RX J1713.7-3946 

2.1. Basic Information from Previous Studies 

The shell-type SNR RX J1713. 7-3946 (G 347.3-0.5), 
which has an angular radius of about 40', was first dis- 
covered during the ROSAT All-Sky Survey (Pfeffermann, 
Aschcnbach 1996). The detection of synchrotron X- 
ray emission with ASCA (Koyama et al. 1997; Slane et 
al. 1999) revealed that the remnant's shell is a presently 
operating source of multi-TcV electrons, which are re- 
sponsible for the synchrotron X-rays. Moreover, the TeV 
7-rays detected by the CANGAROO collaboration from 
this supernova remnant (Muraishi et al. 2000; Enomoto 
et al. 2002) provide direct evidence of the acceleration of 
particles - electrons and/or protons - to energies well be- 
yond 10 TeV. This SNR is very faint at radio wavelengths. 
Based on a lack of thermal X-ray emission, the supernova 
shock is considered to be expanding in a low-density envi- 
ronment, such as a stellar-wind cavity (Slane et al. 1999). 



The distance d, age r, and the explosion energy E§i = 
-E/10 51 erg are three important parameters for under- 
standing the origin of nonthermal radiation components. 
At the perimeter of RX J1713.7— 3946, there is a dense 
(n~lx 10 3 cm~ 3 ) molecular cloud (cloud A) that seems 
to have a very high ratio of two rotational transition 
lines of CO molecules, CO(J=2-1)/CO(J=1-0). This 
suggests an interaction between cloud A and the SNR 
RX J1713.7-3946 (Slane et al. 1999; Butt et al. 2001). 
The distance to cloud A (VLsr= —94 km s _1 ) was kine- 
matically estimated to be 6.3 ±0.4 kpc by using the stan- 
dard Galactic rotation curve (Slane et al. 1999). Thus, if 
the molecular cloud and the supernova remnant are indeed 
physically associated, the distance to RX J1713.7— 3946 
could be significantly larger than the original estimate of 
about 1 kpc based on an absorbing column-density value 
(Koyama et al. 1997). For the preferred distance to the 
source of about 6 kpc, a reasonable set of the age and 
the explosion energy of RX J1713.7— 3946 is in the ranges 
18700-40600 yr and E 51 =2.2-1.7 (Slane et al. 1999). 
However, a different solution of a distance d = 2 kpc, with 
the remnant's age of 2100-3400 yr and the SN explosion 
energy of E51 =1.4-2.2, cannot be ruled out. The smaller 
distance and correspondingly younger remnant agree bet- 
ter with the hypothesis of Wang et al. (1997) who argued 
that RX J1713.7— 3946 could be a likely counterpart of 
the guest star of 393 A.d. in the constellation Scorpius. 

2.2. ASCA Analysis 

The SNR RX J1713.7-3946 was partially observed with 
the ASCA satellite on 1996 September, and most of the 
remnant was covered on 1997 March. The data were 
screened by standard procedures, except that strict rise- 
time screening was not applied to the data taken on 
1997 March. Below we present the results from the Gas 
Imaging Spectrometer (GIS) by adding two identical de- 
tectors (GIS 2 and GIS 3). 

Figure 1 shows the X-ray images of the SNR 
RX J1713.7-3946 in the soft 1-3 kcV and hard 5-10 
keV energy bands. After subtracting the instrumen- 
tal background and correcting vignetting effect, images 
are smoothed with a Gaussian of standard deviation of 
0'.75 and 1 '0 for the soft and hard band, respectively. 
As it has been noted by previous studies (Koyama et 
al. 1997; Slane et al. 1999), shell-like synchrotron X-ray 
images (region 2-4) are seen in both energy bands. In ad- 
dition to region 2-4, the hard-band image clearly shows 
the existence of an extended source (region 1) being co- 
incident with cloud A, although it is barely seen in the 
soft band because of diffuse X-ray emission along the 
Galactic plane. Approximate center of this new source 
is at a celestial coordinate (J2000) of R.A.=17 h 14 m 7?0 
and Decl =-39°12'12"0. We thus refer to this source 
as AX J1714.1-3912. The count rate in the 5-10 keV 
band of AX J1714.1-3912 is 0.029 c s" 1 , which corre- 
sponds to a 34 a statistical significance in excess of the 
background in this region. The total number of source 
photons detected in the 1-10 keV band is derived to be 
3340. Unfortunately, the source lies very near to the edge 
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Table 1. Spectral fits to the ASCA data. 



nonthermal origin. AX J1714. 1—3912 has a good spatial 
association with cloud A (see figure 1). In the following 
discu ssion, therefore, we assume that the X-ray emission 
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1 ' 11 ( 23 Nonthermal particles that are illuminating cloud A in 
72(7^f ^" ra y band could be supplied by either internal or ex- 
98(l4j na ^ acce l era tors. A potential candidate for an external 
72(8S ce l era t° r is the strong shock of SNR RX J1713.7— 3946, 



Note. — Best-fit parameter values and 90% confidence intervals; 
Ajj is equivalent hydrogen column density, Y is photon index, and 
fl_10 keV i s the unabsorbed flux in the 1-10 keV band. \i i s the 
reduced chi-square for the best-fit model and u is the number of 
degrees of freedom of the fit. Region 1-4 are schematically indicated 
by the boxes in figure 1. 

of the observed field. It is therefore difficult to determine 
any extent of the source, which is distorted by the off-axis 
broad point spread function of the telescope. 

The X-ray spectrum of each region (1-4) was derived 
by subtracting the contributions from the Galactic diffuse 
emission taken from a nearby field, and then fitted with 
a photoelectric-absorbed power-law model, described as 
/(e) oc e~ r cxp [— Nntj(e)], where T is a photon index, a(e) 
is the photoelectric cross-section, and iVn is the equivalent 
hydrogen column density. All four spectra were well-fitted 
by these simple models. The best-fit parameters and their 
90% confidence limits arc summarized in table 1. The 
energy spectra and the best-fit models of region 1 (i.e., 
AX J1714. 1—3912) and region 2 are shown in figure 2. The 
power-law indices in region 2-4, T = 2.1-2.3, are consistent 
with previous measurements and are best explained by 
synchrotron radiation from ultrarelativistic electrons. We 
found that the spectrum of AX J1714. 1-3912 is best fitted 
by an extremely flat power-law model with T = LO^'g. 
Within the statistical uncertainties, the absorbing column 
density for region 1 agrees with that obtained for region 
2-4. Attempts at fitting a thermal-bremsstrahlung model 
to the region 1 spectrum gave a lower-limit temperature 
of 35 keV (90% confidence), which is unreasonably high 
for SNRs. 



3. Discussion 

3.1. On the Origin of Flat X-Ray Emission 
We have discovered a new X-ray 
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terized by an extremely flat power-law with T = I.O^q'^. 
The luminosity of the source is estimated to be 
Lx = 1.7 x 10 35 d| erg s _1 , where d 6 is the distance to 
the source normalized to 6 kpc. Since the luminosity 
increases linearly with a high-energy cutoff for the e _1 
type spectrum, the above estimate should be considered 
as a lower limit for the luminosity. It cannot be the 
radiation of thermal gas, unless we assume an extremely 
hot optically thin plasma with a temperature significantly 
exceeding 10 keV. Since the formation of such a thermal 
source seems to be quite problematic, given its extended 
character, below we assume that the X-radiation has a 



which is likely to interact with cloud A, as evidenced 
by a very high ratio of CO(J=2-1)/CO(J=1-0) (Butt 
et al. 2001). The particles can be accelerated in the 
shell of RX J1713.7— 3946 and afterwards enter cloud A, 
given that the shell is a certain site of particle acceler- 
ation, which follows from synchrotron X-ray emissions. 
It is possible that particle acceleration also takes place 
at the secondary shocks inside (or in the vicinity of) 
cloud A initiated by the main shock of the SNR. Finally, 
if there is no physical link between cloud A and SNR 
RX J1713.7— 3946, this would imply the existence of 
unseen internal accelerators inside the molecular cloud. 
Although we cannot exclude any of these possibilities, the 
preferred option seems to be the "cloud interacting with 
SNR" scenario. 

Formally, there are also several options for the non- 
thermal X-ray production mechanisms: synchrotron ra- 
diation, inverse Compton (IC) scattering, and electron 
and/or proton bremsstrahlung. However, here more defi- 
nite conclusions can be drawn. 

3.2. Difficulties in Synchrotron and IC Processes 

Synchrotron X-ray emission by multi-TeV electrons 
could not explain the flat X-ray spectrum. Because of 
the fast synchrotron cooling on a timescale of r syn ~ 
150 (e/5 keV) -1 / 2 yr (for a magnetic field of 25//G typical 
in molecular clouds), the photon index of the synchrotron 
spectrum at X-ray energies becomes T = a c /2+ 1, where 
a c ~ 2-2.2 is the electron acceleration index. Thus, the 
differential spectrum of synchrotron X-radiation cannot 
be flatter than e~ 2 . Note that the power-law fits of the 
X-ray spectra of regions 2-4 (r = 2.1-2.3) agree with the 
synchrotron origin of X-rays from these regions. 

Inverse Compton scattering is precluded by low bright- 
ness of synchrotron radio emission. Relatively low-energy 
(GcV) electrons that produce X-rays via scattering off 
the cosmic microwave background photons, also emit syn- 
chrotron radiation at radio frequencies. Thus, assuming 
that the X-ray flux is due to IC scattering, we can calcu- 
late the radio flux density of the associated synchrotron 
emission. For a magnetic field of 25 fiG the latter is ex- 
pected at a level of 10 4 Jy at 843 MHz, which is higher, by 
three orders of magnitudes, than the upper limit on the 
radio flux from the direction of cloud A (Slane et al. 1999). 
Therefore, we can safely exclude the IC origin of the ob- 
served X-radiation. 

3.3. 1/e Bremsstrahlung 

X-rays can be produced through the bremsstrahlung 
of nonthermal electrons and/or protons interacting with 
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the cloud gas. The energy distribution, for either acceler- 
ated protons or electrons within the cloud, is characterized 
by the position of the break energy, E^, which is set by 
equating the lifetime against the ionization losses and the 
age of the accelerator, Tj on = r agc . At energies below E^, 
for which T lon < r ago , the particle distribution becomes 
"loss-flattened" due to ionization losses (Uchiyama et 
al. 2002). As a result, the bremsstrahlung spectrum obey 
a characteristic energy distribution close to /(e) oc 1/e, 
provided that the energy spectrum of the subrelativistic 
particles becomes harder than E" 1 !" 1 . Almost indepen- 
dent of the details of the acceleration spectrum, this cri- 
terion could be satisfied by the loss-flattened distribution 
below Eb T , which is controlled by the density of the cloud 
n and the source age r age . For n = 10 3 cm -3 and r age = 10 3 
yr, the break appears at E P ^ T ~ 20 MeV in the proton dis- 
tribution, and correspondingly at e Pj b r = (™e/ m p)-Ep,br — 
11 keV in the proton bremsstrahlung spectrum. On the 
other hand, the breaks are E c ^ r = e e ,br — 7 MeV for elec- 
trons. At photon energies anywhere below Ehr, the char- 
acteristic 1/e power-law is predicted. Note that the 1/e 
bremsstrahlung photons below £br are produced predom- 
inantly by particles in a narrow energy interval around 
i?br- The spectral shape of the flat X-ray spectrum of 
AX J1714. 1—3912 agrees perfectly with the 1/e emis- 
sion of either proton bremsstrahlung (PB) or electron 
bremsstrahlung (EB) from the loss-flattened distribution. 

Whereas the PB and EB mechanisms are expected to 
give rise to the same 1/e X-ray spectrum, the PB luminos- 
ity peaks at (hard) X-rays and the EB peaks at gamma- 
rays for typical parameters. In the case of the 1/e EB 
emission, the photon flux per each logarithmic interval 
in photon energy anywhere up to e e ,br ~ 7 MeV is con- 
stant. Then, the observed X-ray flux of 4 x 10 -3 photon 
cm~ 2 s _1 in the 1-10 keV band implies that the 1-10 MeV 
flux should be comparable to the Crab. This is incon- 
sistent with non-detection by the COMPTEL instrument 
(Schonfelder et al. 1996) onboard the Compton Gamma 
Ray Observatory. Moreover, unless we introduce a sharp 
cutoff in the electron distribution at an energy of around 
100 MeV, the EB model overshoots the 7-ray flux reported 
by the EGRET instrument. Therefore, we may conclude 
that the gamma-ray data favor proton bremsstrahlung, 
rather than the electron bremsstrahlung scenario. 

3. 4- Energetics Requirement of 1/e Bremsstrahlung 

The X-ray luminosity of the 1/e bremsstrahlung emis- 
sion of either protons or electrons can be estimated as 

Lx ~ ?7p,o — — , (1) 

\ 7"agc / 

where W is the total amount of kinetic energy contained 
in X-ray emitting particles, T agc is the age of the source 
(accelerator), and 77 = Tbn/Tbrem- Here, we take into ac- 
count that the bulk of X-rays are produced by particles 
from the break region, thus T agc = T lon . For the parameters 
chosen above, ?y p = 3 x 10 -5 and ry = 9 x 10 -5 . These fac- 
tors demonstrate the low efficiency of the bremsstrahlung 
mechanism relative to the ionization losses, implying a 



huge, W p = 5.6 x 10 39 d| erg s" 1 and W c = 1.9 x lO 39 ^ 
erg s _1 , injection power in protons and electrons, respec- 
tively Correspondingly, the total amounts of energy re- 
leased during operation of the accelerator are Wp Uh = 1.8 x 
10 50 d|(r ago /10 3 yr) erg and W c = 7.2 x 10 49 ^( Tagc /10 3 yr) 
erg. 

If the particles are accelerated in the shell of SNR 
RX J1713. 7-3946, and only a relatively small (10% or so) 
fraction of these particles enter cloud A, the total kinetic 
energy, which is transferred to subrelativistic protons, 
should be at least 10 51 erg; the electron bremsstrahlung 
model alleviates the total energy by a factor of about 3. 
Given the limited energy budget of a supernova shock of 
about 10 51 erg, we face a serious problem to support the 
required X-ray luminosity, unless the system is closer to 
the Earth than estimated by the radial velocity of cloud A. 
For example, a distance of 2 kpc would reduce the energy 
requirement to a quite comfortable level of ~ 2 x 10 50 erg. 
Several other scenarios may be invoked to overcome the 
energy budget problem. The simplest assumption would 
be that the explosion energy significantly exceeds 10 51 
erg, and that more than 10% of this energy is released 
in the form of subrelativistic particles. It is also possible 
that the nonthermal particles are accelerated inside the 
cloud, e.g. by the shock initiated at the collision of the 
blast wave of RX J1713. 7-3946 with cloud A (Bykov et 
al. 2000). Finally, we may speculate that the essential part 
of the energy which goes into the ionization and heating of 
gas, and also the excitation of plasma waves, is returned 
to subrelativistic protons through acceleration on these 
plasma waves. Apparently, all of these assumptions need 
detailed quantitative studies. 

AX J1714. 1-3912 and cloud A are located within 
the error box of the unidentified 7-ray source 3EG 
J1714-3857 (Hartman et al. 1999). The GcV 7-ray flux 
may be explained by 7r° gamma-rays by collisions be- 
tween the shock-accelerated protons and cloud A. Butt 
et al. (2001) argued that an electron-bremsstrahlung ori- 
gin for the GeV flux is less likely because of the faint 
synchrotron radio emission of cloud A. The total en- 
ergy liberated in 7-rays through the decay of neutral 
pions is roughly estimated to be L 1 ~ (l/3)o- pp ncWp Cl 
where cr pp ~ 30 mb is the proton-proton inelastic cross- 
section, and W p el is the energy content in GeV protons, 
within cloud A. The EGRET luminosity, L 1 = 6.5 x 10 35 d| 
erg s" 1 , corresponds to W*° l = 2.2 x 10 48 d|(7i/10 3 c m - 3 )- 1 
erg. Thus, with the proton-bremsstrahlung model for 
AX J1714. 1—3912, the energy content in subrelativistic 
protons far exceeds that in relativistic protons, W p uh ~ 
80 W p cl . In this case, it would appear that the accelera- 
tion/injection mechanisms allow a small fraction of pro- 
tons to be accelerated to relativistic energies. Another 
possibility could be that the bulk of accelerated relativistic 
protons have already diffused away, while the subrelativis- 
tic protons are still being captured within the molecular 
cloud. 
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Fig. 1. X-ray surface brightness (color map linearly spaced 
from black to white) of the SNR RX J1713.7-3946 in (a) 
1-3 keV and (b) 5-10 keV energy bands obtained with the 
ASCA GIS detectors. Right ascension and declination are 
in the J2000 coordinate system. The white thick line repre- 
sents the Galactic plane. A transient X-ray source (Koyama 
et al. 1997) is blanked out [a white cross in (b)]. Depicted are 
the 68% confidence error contour of the EGRET unidentified 
source 3EG J1714-3857 (Butt et al. 2001), and an approx- 
imate outline of cloud A based on the CO(J=l— 0) intensity 
(Bronfman et al. 1989). 
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Fig. 2. Unfolded X-ray spectra of region 1 
(AX J1714.1-3912) and region 2 measured by the ASCA 
GIS detectors. The solid lines are the best-fit power-law 
models, taking account of the photoelectric absorption at 
lower energies (see table 1), while the dashed lines are 
absorption-corrected power-law spectra. For illustrative 
purposes, we show the X-ray spectra of region 1 being derived 
by using an alternative background field (open circles). 



